Las Pailas geothermal field - Central America case study: Deciphering a volcanic geothermal play type through the combination of optimized geophysical exploration methods and classic geological conceptual models of volcano-tectonic systems

Sustainable exploitation strategies of high-enthalpy geothermal reservoirs in a volcanic geothermal play type require an accurate understanding of key geological structures such as faults, cap rock and caldera boundaries. Of same importance is the recognition of possible magmatic body intrusions and their morphology, whether they are tabular like dikes, layered like sills or domes. The relative value of those magmatic bodies, their age, shape and location rely on the role they play as possible local heat sources, hydraulic barriers between reservoir compartments, and their far-reaching effect on the geochemistry and dynamics of fluids. Obtaining detailed knowledge and a more complete understanding at the early stages of exploration through integrated geological, geophysical and geochemical methods is essential to determine promising geothermal drilling targets for optimized production/re-injection schemes and for the development of adequate exploitation programs. Valuable, extensive geophysical data gathered at Las Pailas high-enthalpy geothermal field at northwestern Costa Rica combined with detailed understanding of the geological structures in the underground may represent a sound basis for an in-depth geoscientific discussion on this topic. Currently, the German cooperation for the identification of geothermal resources in Central America, implemented by the Federal Institute for Geosciences and Natural Resources (BGR), supports an international and interdisciplinary effort, driven by the Instituto Costarricense de Electricidad (ICE) with different international and national research institutions, including the Leibniz Institute for Applied Geophysics (LIAG). The discussions and joint studies refer to the optimized utilization of geophysical and geological methods for geothermal exploration in the Central American region, using the example of Las Pailas Geothermal Field. The results should contribute to a better understanding of the most appropriate geothermal exploration concepts for complex volcanic field settings in Central America

Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and investigate numerically the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geologocal configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems

Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and investigate numerically the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geologocal configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems

Bilinear pressure diffusion and termination of bilinear flow in a vertically fractured well injecting at constant pressure

This work studies intensively the flow in fractures with finite hydraulic conductivity intersected by a well injecting or producing at constant pressure, either during an injection or production well test or the operation of a production well. Previous investigations showed that for a certain time the reciprocal of flow rate is proportional to the fourth root of time, which is characteristic of the flow regime known as bilinear flow. Using a 2D numerical model, we demonstrated that during the bilinear flow regime the transient propagation of isobars along the fracture is proportional to the fourth root of time. Moreover, we present relations to calculate the termination time of bilinear flow under constant injection or production well pressure as well as an expression for the bilinear hydraulic diffusivity of fractures with finite hydraulic conductivity. To determine the termination of bilinear flow regime, two different methods were used: (a) numerically measuring the transient flow rate in the well and (b) analyzing the propagation of isobars along the fracture. Numerical results show that for low dimensionless fracture conductivities the transition from bilinear flow to another flow regime (e.g., pseudo-radial flow) occurs before the pressure front reaches the fracture tip, and for high dimensionless fracture conductivities it occurs when the pressure front arrives at the fracture tip. Hence, this work complements and advances previous research on the interpretation and evaluation of well test analysis under different reservoir conditions. Our results aim to improve the understanding of the hydraulic diffusion in fractured geologic media, and as a result they can be utilized for the interpretation of hydraulic tests, for example to estimate the fracture length

Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and numerically investigate the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and a variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geological configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems.Leibniz-Gemeinschaft https://doi.org/10.13039/50110000166

Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and numerically investigate the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and a variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geological configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems

Pressure transient analysis to investigate a coupled fracture corridor and a fault damage zone causing an early thermal breakthrough in the North Alpine Foreland Basin

The heterogeneity of the Upper Jurassic carbonate reservoir (Malm reservoir) beneath the North Alpine Foreland Basin has a significant influence on the mass and heat flow processes during geothermal exploitation. Geophysical borehole data revealed that sub-seismic scale fractures and karstified fractures occur at the inflow zones of deep geothermal wells. However, pressure transient analysis (PTA) in some previous studies concluded that it is difficult to detect the influence of sub-seismic scale features, suggesting that radial flow regime is dominant. Accordingly, a regional thermal-hydraulic model adopted the equivalent porous medium (EPM) approach, homogenizing the sub-seismic scale reservoir heterogeneities; however, unable to detect an early thermal breakthrough (ETB) in a geothermal doublet located SE of Munich. We apply PTA on three buildup tests belonging to that doublet following a deterministic approach to constrain the reservoir type by interpreting the pressure derivative (PD) plots constrained by geophysical and geological data. We derive the magnitudes of the reservoir hydraulic parameters by matching the PD plots with the selected interpretation models. We find that clustered fractures have a significant influence on the reservoir hydraulics, evidenced by trough-shaped curves in the PD plots. Linear flow regime interpreted from the interference test between the two wells indicates permeability anisotropy, which may have caused the ETB. Geophysical data interpretations indicate that these fractures correspond to a coupled fault damage zone and a fracture corridor. Finally, we present a fit-for-purpose 2D discrete fracture network model utilizing the PTA results to match our analytically calibrated model. Our study offers a potential hydraulic explanation to the cause of the ETB highlighting the importance of integrating multi-scale/disciplinary data sets to improve the reliability of dynamic reservoir models, based on which, economic-related decisions are made.Applied Geolog